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Abstract 

It is shown experimentally that vertical pairing of two identical microspheres 
suspended in the sheath of a radio-frequency (rf) discharge at low gas pres- 
sures (a few Pa), appears at a well defined instability threshold of the rf power. 
The transition is reversible, but with significant hysteresis on the second stage. 
A simple model, which uses measured microsphere resonance frequencies and 
takes into account besides Coulomb interaction between negatively charged 
microspheres also their interaction with positive ion wake charges, seems to 

explain the instability threshold quite well. 
PACS number(s): 52.25. Zb, 52.25.Ub 



1 



Typeset using REVTgX 



Recent interest in the properties of complex plasmas - charged monodisperse micropar- 
ticles suspended in an electron-ion environment, is partially due to the possibility to model 
condensed matter phenomena on an "atomic" level The particles are charged nega- 

tively in a radio-frequency (rf) discharge (up to 10^ — 10^ elementary charges) and levitate 
in the lower plasma sheath, where gravity can be balanced by an inhomogeneous vertical 
electric field. The particles repel each other via a screened Coulomb interaction. Inside a 
confining horizontal potential they arrange themselves either into an ordered structure - the 
"plasma crystal" , or form disordered, liquid- or gas-like, states . 

At low gas pressures one has typically only a very few layers (as few as one) for these 
structures. Ion streaming motion in the plasma sheath produces a nonequilibrium environ- 
ment, which is reflected in the properties of the complex plasma structures and makes them 
rather different from the known equilibrium ones One of the striking observations at 

sufficiently low gas densities (such that the ion-neutral mean free path is not small compared 
with a typical distance between microspheres) is an unusual "stacking" of the particles such 
that adjacent horizontal layers are located on top of each other 0,^,^. This vertical "po- 
larization" of the dusty plasma was ascribed to the attractive forces which result from the 
focusing of the ion flux under each particle (ion wake effect, see, e.g., Ref. []TD[)- Thus, the 



attractive interaction between the horizontal layers in a plasma crystal is asymmetric, such 
that the attractive force is communicated only downstream in the direction of the ion flow 
T0| . The same arguments were applied to predict a "binding" force for a possible "dust 



molecule" formation |TT[. This theoretical prediction was verified recently by experimental 
studies of the competition between repulsive and attractive forces acting on two different 
dust particles levitated on different levels in the plasma sheath [jl2 . 



In this Letter we present experiments on a new instability which is observed for two 
identical microspheres suspended initially on the same level in the plasma sheath. In this 
work we will report on the quantitative investigations involving (for simplicity) only two 
microspheres, which were undertaken to clarify the basic behavior for many identical par- 
ticles. The instability appears first as a continuous (forward) bifurcation to a state of the 



vertically separated microspheres, and then to a discontinuous vertical pairing. Depending 
on gas pressure, the last stage - the pairing - can be either strongly hysteretic (for lower 
pressure) or weakly hysteretic (for higher pressure). 

The experiments were performed in a standard Gaseous Electronics Conference (GEC) rf 
Reference Cell |TB| with the lower electrode powered at 13.56 MHz and the upper electrode 
grounded (see Fig. p. Argon gas at various pressures between 1 and 7 Pa was used, and 
the discharge power (or a rf peak-to-peak voltage, Upp) was the control parameter of the 
instability described. The electron temperature and density were measured at the center 
of the discharge with a rf-compensated passive Langmuir probe. At these conditions the 
electron temperature was found to be within 2-5 eV, while the electron density Ug ranged 
from 10'' to 8 X 10*^ cm~^, so that the electron Debye length varied from Ade — 0.5 mm (high 
Upp) to — 5 mm (low Upp). These values of the Debye length are comparable to, or larger 
than the separation distances measured, and thus Debye shielding will play only a minor 
role in this rather low-density plasma. The microspheres, suspended in the plasma, were 
illuminated by a laser sheet of about 100 /im thickness, and their imaging was performed by 
external CCD cameras from the top and from the side (in both cases via 300 mm objectives). 
We used polystyrene (density 1.05 g/cm^) microspheres with diameter 7.6±0.1 fim and mass 
M ~ 2.4 X 10^^° g. In order to confine the microspheres horizontally we placed a quartz 
glass cylinder of 20 mm height and 50 mm diameter on the lower electrode (see Fig. p. 
In contrast to a copper ring normally used to confine particles, the glass cylinder produces 
extremely flat, almost "square well" confining potential. The latter conclusion is based on 
our observation of unusually fiat shape of 2D plasma crystal, particularly up to its edges, 
contrary to the trap formed by a standard metal ring. The "square well" confining potential 
provides us the opportunity to study two identical microspheres suspended in the plasma 
sheath on the same initial height. 

In order to quantitatively investigate the pairing instability and to verify that it is indeed 
symmetry breaking, we conducted experiments with two identical microspheres at different 
pressures. Figure ^ shows, how, at a pressure of 2 Pa, the relative position of one microsphere, 
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which eventually moves below and then under the other (i.e., downstream) with respect to 
the eventual upstream microsphere, varies as the power (or f/pp) decreases. As long as Upp 
exceeds the threshold value ~ 55 V, both particles are located on the same level (initial 
position 1). When Upp decreases below the threshold, vertical separation starts (position 
2) and grows continuously up to position 3. A further small Upp decrease causes the lower 
particle to "jump" to a position directly beneath the upper one and to create a vertical pair 
(final position 4). By reversing the process, i.e., increasing Upp, the particles revert back to 
their initial horizontal configuration. The transition is strongly hysteretic with respect to 
Upp. Further measurements showed that the threshold, U^, rapidly increases with pressure, 
while the relative width of the hysteresis decreases with pressure. 

Figure ^represents the forward (positions 1-5) and the reverse (positions 6,7) transitions 
aX p = 7 Pa. The instability (vertical separation) starts at U^ ~ 200 V and develops 
continuously until position 4. Then a small Upp decrease is accompanied by the discontinuous 
vertical pairing (position 5). When the voltage is increased back, the particles remain 
vertically paired until position 6 is reached, and then the lower particle "jumps" to position 7. 
In Fig. ^ the vertical separation distance is plotted against the control parameter, U^ — Upp, 
at p = 7 Pa. As in the low pressure case, the transition is divided into two stages: (i) 
continuous transition from the horizontal configuration, and (ii) discontinuous, hysteretic 
transition to the final vertical pairing. It is worth noting that at the onset of the continuous 
transition strong symmetrical fiuctuations of the vertical particle separation are observed 
in the experiments. This means that the pairing instability may be initiated with equal 
probability by either microsphere in the pair. Therefore, this transition is related to the 
symmetry breaking. 

In the experiments both Upp (at the fixed p) and p (at the fixed Upp) were used as 
the control parameters of the instability. But in a subsequent analysis, presented below, it 
is much more convenient to use the resonance frequencies of vertical, ujz, and horizontal, 
ujr, oscillations of a single particle in the sheath as the control parameters. The resonance 
frequencies are certain functions of Upp and p, and can be determined experimentally . 



We first examine the simplest model without wake effects. The model uses the measured 
values of uj^ and as the control parameters and provides an order of magnitude estimate 
for the instability threshold. We consider a pair of two identical particles of mass M and 
(negative) charge — Q, electrostatically confined in the harmonic potential well. If the par- 
ticles are separated horizontally by R and vertically by 5 (i.e., each particle is displaced at 
R/2 and 5/2 from the center), then their energy in the confinement is ^Miuld"^ + uf.R'^). 
Since the Debye length is relatively large, the bare coupling energy of two unshielded par- 
ticles, Q^{R^ + 5^)"^''^, should be added to the confinement energy. This system has two 
stable configurations - horizontally (when > oOr) or vertically (when uj^ < ujr) aligned at 
distances of R = {2Q'^ / Mulfl'^ oi 5 = {2Q'^ / Mu'j^)^'^ , respectively. 

Our model is characterized in terms of the measured resonance frequencies. Thus, it 
is necessary to trace these frequencies as the control parameter, U'^ — f/pp, is varied. The 
typical dependence of Uj. and Uz on the value of f/pp at p = 1 Pa is shown in Fig. |^. 
One can see that the frequencies converge rapidly as f/pp decreases. Comparison of Figs. ^ 
and 1^ shows (albeit for different pressures) that the transition starts continuously when 
ujz still considerably exceeds uo^. This implies that the harmonic approximation with the 
isotropic particle field is too simple. Thus, we have to include the wake effect in its simplest 
form. For consistency in the nonlinearity analysis, we must also take into consideration 



the anharmonicity of the potential well in the vertical direction |]T5[. The corresponding 



expression for the total potential energy of a particle pair in a 2D confinement is 

Wpa. ^ - i U^R' + 0015' + ^5'\ + 



4 V ' ' 8 / V^TP 2^W+{A + 6Y 2^R^ + {A-6)^' 

(1) 

In writing Eq. (|I]) the simplest two- microsphere model of the wake is used (see Fig. ^). We 
treat the excessive positive charge of the wake, q, as point-like, located under the particle 
at distance A. We assume also that since both the particle and the wake charges vary only 
slightly with variation in R and 6, so that Q and q are approximately constant. The last two 
terms in Eq. (|l|) represent the particle-wake interaction, and cross-terms of the particle-wake 



and wake- wake (gg) interactions are neglected [|T2[. Since the coefficients multiplying the 



terms 6^ cancel out for a pair of identical particles, the lowest order vertical anharmonicity 
in Eq. (|1|) is the fourth order term with coefficient (3* > but whose value is here unknown, 



since it was not feasible to measure it by technique used elsewhere [0. 

The dependence R{S) is determined from the equilibrium condition in the radial direction, 

9Wpair/9i? = 0: 

2Q^ / 1 1 \ 2 

(i?2 + ^2)3/2 - Ql (^[^2 + (A + 5)2]3/2 + [^2 + (A - 5)2]3/2 j = ^"^r- (2) 

The equilibrium value of 5 is given by the condition dWpair/d^ = 0. A linear combination 
of this condition and Eq. (Q) gives 

MB* / 1 1 \ 

- -.^)« + — «^ + i^^j-^^-^ - j = 0. (3) 

Expanding Eq.(^ to third power in 5 and substituting R{5) from Eq. (Q) provide a typical 
stationary equation for the vertical displacement as the order parameter: 

V"^-' V 6(1 + A'^)^ ^ 1 - ) 

Here we have introduced the dimensionless parameters, 5 = 5/Rq and A = A/Rq [normalized 
by -Ro = R\6=o from Eq. (H)], and the critical frequency ratio. 



Q = 



^ 6ggA2 



The instability criterion {uJz/uJr)cT = is obtained from Eq. by setting the coefficient 
of 6 equal to zero. At Uz/ur > fi, the energy Wpair(5) has a minimum at 5 = and the 
equilibrium configuration is horizontal. In the opposite case, oJz/uOr < the 5 = state 
becomes unstable, and the particles start separating vertically. From Fig. ^ we note that 
the separation starts at {ijJz/^jJr)cT — 4 (when t/pp ^ f^pp)- Assuming Q ^ q and substituting 
Eq. (0) at 5 = into Eq. (^, one gets f2 as a function of A only. Then at A ~ 1 we obtain 
f2 ~ 3. Thus, the threshold of the pairing instability is adequately described by the model. 
A saturation of the order parameter, 5, above the instability threshold in the ffist continuous 
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stage of the transition (see Fig. ^ could be reached if the sum of both coefficients of in 
Eq. @ is positive. By using the expression f2(A) and the value of the ratio Q/g ~ 1 we 
achieve this condition at (3*R^/uj'j ^ 7, or /3*/(27r)^ ^ 20 Hz^/mm^. Our recent experiments 



15| indicate that at pressures below ~ 10 Pa the anharmonic coefficients become relatively 



large: For example, /5/(27r)^ ~ 20 Hz^/mm^ at p ~ 1 Pa and f/pp ~ 70 V. At smaller f/pp, 
the whole vertical structure of the sheath changes dramatically, because the Debye length 
increases by the order of magnitude. Hence, the anharmonic coefficients should also increase 
rapidly as U^^ decreases, and it is plausible that (3* will provide the instability saturation. 



Then, as follows from Eq.(|D, 6 oc — oj^/^r close to the threshold. However, failing 
an accurate experimental determination of the shape of the sheath potential for the actual 
conditions of the experiment, further testing of the model via the experimental value of 5 is 
not feasible here. 

The modelling of the second (discontinuous) stage of the instability is rather complicated 
problem: When the lower particle approaches the wake, its actual "shape" becomes crucial, 
and the approximation of the point-like wake charge cannot be used. This very nonlinear 
second stage will almost certainly involve nonlinear wake dynamics, and these can probably 
only be adequately treated by three-dimensional numerical simulations. We would like to 
point out that the effect of the wake is stronger at lower p and t/pp, when Adb is larger. 
At the same time, the ion mean free path should exceed Adb (otherwise, the wake effect 
weakens due to the ion scattering on neutrals). Both these requirements are well satisfied 
in our experiments. 

V. S. is grateful for the support of the Alexander von Humboldt Foundation, Germany, 
during his stay at MPE. Many useful critical remarks of an anonymous Referee are gratefully 
acknowledged. 
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FIGURES 



FIG. 1. Experimental setup using a quartz glass cylinder as a confining electrode. 

FIG. 2. Sequence of the relative particle positions during the pairing instability as Upp, 
decreases (pressure p = 2 Pa). The steps are: (1) Upp > ~ 55 V - horizontal configuration, (2) 
Upp ~ 45 V and (3) Upp ~ 37 V - continuous separation, (4) Upp ~ 36 V - discontinuous vertical 
pairing. 

FIG. 3. Sequence of the relative particle positions during the pairing instability as Upp, de- 
creases (1-5) and increases back (6,7) (pressure p = 7 Pa). For decreasing voltage the steps are: 
(1) Upp > U^l ~ 200 V - horizontal configuration, (2) Upp ~ 110 V, (3) Upp ~ 85 V, and (4) 
Upp ~ 78 V - continuous separation, (5) Upp ~ 76 V - discontinuous vertical pairing; For increas- 
ing voltage the steps are: (6) Upp 2± 93 V - continuous vertical movement and (7) Upp ~ 95 V - 
discontinuous reversibility of the pair configuration. 

FIG. 4. Vertical separation of a particle pair vs. the control parameter, Uj^ — Upp, for two 
stages of the pairing transition at p = 7 Pa. The hysteretic nature of the second stage is obvious 
from the plot. The arrows denote the direction of the variation observed. 

FIG. 5. Horizontal, ujr {o), and vertical, ujz (+), resonance frequencies of a single particle as 
function of Upp (at a pressure p = 1 Pa). The vertical dashed line shows the threshold voltage 
C/pp ~ 35 V (onset of the vertical separation) . 

FIG. 6. Schematics of the simplified model of the wake potential. 
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